Localization of Ultrasound in a Three-Dimensional Elastic Network* 
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After exactly half a century of Anderson 
localization 1 , the subject is more alive than ever. 
Direct observation of Anderson localization of 
electrons was always hampered by interactions 
and finite temperatures. Yet, many theoret- 
ical breakthroughs were made, highlighted by 
finite-size scaling 2 , the self-consistent theory^ and 
the numerical solution of the Anderson tight- 
binding modePES. Theoretical understanding is 
based on simplified models or approximations 
and comparison with experiment is crucial. De- 
spite a wealth of n ew experi mental data, with 
microwave s^^, lig ht)^- 9 * 10 * 11 * 12 ^ ultrasound^ and 
cold atoms^ 1 * 15 * 16 *, many questions remain, es- 
pecially for three dimensions. Here we report 
the first observation of sound localization in a 
random three-dimensional elastic network. We 
study the time-dependent transmission below the 
mobility edge, and report "transverse localiza- 
tion" in three dimensions, which has never been 
observed previously with any wave. The data 
are well described by the self-consistent theory 
of localization. The transmission reveals non- 
Gaussian statistics, consistent with theoretical 
predictions. 

Most text books on condensed matter explain that the 
electronic states in disordered conductors are extended 
plane or Bloch waves with finite life times. This gives rise 
to "Ohmic resistance" , proportional to the length of the 
sample. In the picture presented by AndersorP, "large" 
disorder makes electronic states localized in space. This 
offers a mechanism to explain the widely observed metal- 
insulator transitions 17 . Scaling theory proposes a sin- 
gle parameter, the Thouless conductance <?, to describe 
the anomalous length dependence of the resistance of a 
sample^. The localized regime is characterized by the so- 
called Thouless criterion g < 1. While these ideas were 
first proposed for electron localization, in the early eight- 
ies interest in classical wave localization was raised^ES 
with the promise of avoiding the difficulties caused by 
interactions in electronic systems. At the same time, the 
absence of bound states for classical waves makes local- 
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FIG. 1: The random elastic network and correspond- 
ing transmission coefficient, a, Top view of one of the 
samples (right), along with a magnified image showing the 
mesoscale structure of the random elastic network (left). The 
samples were made from 4.11 ± 0.03 mm diameter aluminum 
beads brazed together at a volume fraction of approximately 
55%. The diameter (120 mm) of the slab-shaped samples was 
much larger than the thickness L, which ranged from 8.3 mm 
to 23.5 mm. b, Frequency dependence of the amplitude trans- 
mission coefficient for L = 14.5 mm. The arrows indicate the 
resonance frequencies of isolated aluminum beads. 

ization more challenging to achieve in practice, with ab- 
sorption as an additional concern. 

In this Letter we demonstrate Anderson localization 
of ultrasound in a three-dimensional (3D) medium. Our 
samples are single-component random networks made by 
brazing aluminum beads together, see Fig. [TJl. With 
ultrasound we probe the vibrational excitations of the 
network in the intermediate frequency regime (0.2 to 3 
MHz), where the wavelength is comparable to the bead 
and pore sizes. We use pulsed techniques to measure the 
amplitude transmission coefficient, shown in Fig. [TJj. The 
transmission spectrum exhibits a succession of band gaps 
and pass bands, due to the overlapping resonances of the 
aluminium beadiP^. Our study focuses on frequencies 
just below the first band gap at 0.5 MHz, and at higher 
frequencies (1.6-3 MHz), where it was possible to extract 
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the phase and amplitude of the coherent pulse for longi- 
tudinal waves crossing the sample. Hence we were able 
to measure the longitudinal phase v p and group v a ye- 
locities, as well as the scattering mean free path pHHl 
Note that since absorption is negligible (see below), the 
attenuation of the coherent pulse gives £ directly. Around 
0.20 MHz, v p = 1.75 km/s, v g = 2.1 km/s and £ = 2.2 
mm, while around 2.4 MHz, v p — 5.0 km/s, v g = 5.2 
km/s and £ = 0.6 mm. This leads to a product of wave 
vector k and scattering mean free path k£ = 1.6 at 0.20 
MHz and k£ — 1.8 at 2.4 MHz. The values of fc^ for shear 
waves, likely to dominate, are not known but are proba- 
bly roughly equals The small values of kl indicate that 
our samples are strongly scattering. For 3D disorder lo- 
calization is expected when kl < l^. Because the exact 
critical value is not known, the ultrasonic waves at both 
frequencies are candidates to be Anderson localized. 

In previous reports on Anderson localization with clas- 
sical waves, absorption has been a major obstacle to 
reaching unambiguous conclusion d 7 * 8 * 10 * 11 * 25 !. The fol- 
lowing experiment is capable of probing Anderson local- 
ization without being blurred by absorption. We mea- 
sure the spatially and time-resolved transmitted inten- 
sity through our sample. Using a quasi-point source 
that is about a wavelength wide and a sub-wavelength- 
diameter detector that scans the sample at various trans- 
verse positions p in the near field (p = opposite to 
the source) , the average transmitted intensity I(p, t) was 
determined^. From these measurements, we determine 
the ratio I(p, t)/I(0,t), probing the dynamic spreading 
of the intensity in the transverse direction. Any pos- 
sible absorption factor exp(— t/r a ) cancels in this ra- 
tio. For each p, a transverse width w p (t) of I(p,t) can 
be defined by setting I(p, t)/I(0, t) = exp[— p 2 / w 2 p (t)]. 
If the wave propagation is diffuse, the spatial intensity 
profile is Gaussian and w p (t) = 4Dt is independent of 
p. We have observed this normal diffuse behavior in 
less strongly scattering samples, and hence been able to 
demonstrate a way of measuring the diffusion coefficient 
without the usual complications due to boundary condi- 
tions and absorption^. 

For frequencies between 1.6 and 3 MHz, we observe 
dramatically different behavior, shown in Fig. [2] for two 
representative samples at 2.4 MHz. Instead of increasing 
linearly with propagation time, w 2 {t) is seen to saturate. 
This saturation is reminiscent of transverse localization, 
previously predicted in 3D systems with 2D disorder-^ 
and observed in 2D disordered photonic crystals^. In 
our samples, the disorder is clearly 3D since the scat- 
tering mean free path is much smaller than the sample 
thicknesses (10 < L/£ < 40). It is not at all clear a 
priori that transverse localization can occur with 3D 
disorder. Scaling theory^ predicts a diffusivity D{L) de- 
pendent on the sample size L, so that w 2 (t) — 4D(L)t 
still rises linearly with time. The saturation of w p (t) 
could possibl y be ex plained by a diffusivity D(t) decay- 
ing with tirn d 9 ! 10 ^ . This would still lead to a Gaussian 
transverse intensity profile and hence to w 2 constant with 
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FIG. 2: Transverse localization in three dimensions. 

a,b, Temporal evolution of the effective width squared, Wp(t), 
of the transmitted intensity emanating from a point source for 
several transverse displacements p of the detector and for two 
sample thicknesses. The frequency is 2.4 MHz. The solid 
curves are the best fits of our self-consistent theory to the 
experimental data (symbols), from which we determine i* B = 
2 mm, Db = 11 m 2 /s, £ « 15 mm for the L = 14.5 mm sample 
and £ « 7 mm for the L = 23.5 mm sample. Other input 
parameters — kl = 1.8 and the internal reflection coefficient 
R = 0.82 — were obtained from independent measurements. 
The dashed line shows the linear time-dependence of w 2 p that 
would be seen for diffuse waves, using D = 1.25 m 2 /s. c, 
Dependence of the intensity ratio on distance p at six different 
times, showing the non-Gaussian profile that is found both 
experimentally (symbols) and theoretically (solid curves). 



p. However, this is not what we observe. 

To describe the dynamics of the anomalous transverse 
confinement of the multiply scattered waves, we apply 
a novel version of the self-consistent (SC) theory of lo- 
calization. The new element consis ts in incorporating 
boundary conditions self-consistently^^H. This theory 
provides a position dependent, dynamic diffusivity ker- 
nel D(r,t — t'). Near the mobility edge, the position 
dependence of D affects all aspects of wave transport 
considerably. The SC theory requires as input the value 
for k£, the localization length £, the diffusion constant 
Db free from macroscopic interferences, and the internal 
reflection coefficient R. In the model we replace the inci- 
dent focused beam by a point source at depth £* B . This 
is the transport mean free path associated with diffusion 
in the absence of macroscopic interferences, which ought 
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to be negligible just after the beam comes in. In Fig. [2] 
we compare the observed dynamics of transverse width 
to this theory. Excellent agreement with the data is seen 
for all p with a single set of parameters for each sample 
(solid curves), yielding, in particular, £ ss 15 mm and 
£ 7 mm for the thinner and thicker samples respec- 
tively. We attribute the smaller value of £ for the thicker 
sample to stronger scattering due to small differences in 
the microstructure. This is consistent with the strong 
sensitivity of £ to small changes in disorder that is ex- 
pected near the localization transition. As p increases, 
the curves w 2 p {t) move upwards, meaning that the ob- 
served transverse profile I(p,t) is not Gaussian. Figure[2] 
shows that this behavior is well captured by the theory, in 
which the position dependence of D is a crucial element. 
Any homogeneous absorption would not affect Fig. [2] We 
believe that this combination of theory and experiment 
provides strong evidence for Anderson localization of ul- 
trasound near 2.4 MHz in our samples and allows us to 
estimate the mobility edge (k£) c , which we find to be 
approximately 1% above the measured k£ = 1.8. 

To find additional support for our conclusions, we have 
measured the time-dependent transmission using an ex- 
tended quasi plane-wave source. Near 0.2 MHz, the 
average transmitted intensity I(t) was found to decay 
exponentially at long times (Fig. [3^), with the entire 
time dependence of I(t) being well described by diffu- 
sion theory^. We conclude that multiply scattered ul- 
trasound at 0.2 MHz propagates in a normal, diffuse way. 
By contrast, for ultrasound propagating near 2.4 MHz, 
the time dependence of I(t) shows a quite different be- 
haviour (Fig. ^ijp), with a non-exponential tail at long 
times. Similar behaviour has been reported before and 
is often explained by a time-dep endent reduction of the 
effective diffusion coefficient £)(^ 9 I 10 I11J The solid curve 
in Fig. [3]} is a fit to the SC theory, and gives an ex- 
cellent description of the experiment at all propagation 
times. The good agreement between theory and experi- 
ment supports our previous conclusion that ultrasound at 
2.4 MHz is localized. From the fits in Figs. [2]and|3]3 and 
the relation Db — ve£*b/3 w ^ find transport velocities 3 
to 5 times larger than the phase velocity of longitudinal 
waves. Further theoretical work is needed to understand 
these — apparently large values for ve in the localized 
regime. 

We also address the statistical approach to Ander- 
son localization^. The normalized transmitted intensity 
1 = 1/ (I) was measured for a large number of individual 
speckle spots in the near field, and for a broad incident 
beam. The intensity distributions P(I) for both frequen- 
cies are shown in Fig. [4] We have compared our data to 
the theory of Ref!™, with the Thouless conductance g as 
the only free parameter in the fit. The agreement is excel- 
lent for g = 0.80 ± 0.08 at 2.4 MHz and g = 11.4 ± 0.8 at 
0.2 MHz. The strong deviation from Rayleigh statistics 
with g < 1 observed at 2.4 MHz is interpreted as a sig- 
nature of localization^. The remarkable agreement with 
theory, derived for the intensity in the far field, and for 
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FIG. 3: Averaged time-dependent transmitted inten- 
sity. Transmitted intensity, I(t), normalized so that the peak 
of the input pulse is unity and centered on t — 0, at represen- 
tative frequencies in the diffuse (a) and localized (b) regimes 
for the L = 14.5 mm sample investigated in Figs. 2a, c. In a, 
the fit to the diffusion theory with R — 0.85 gives D = 3.0 
m 2 /s, I* ~ 2.5 mm and r a too large to be measurable. In b, 
the data are fitted by the SC theory (red curve), with £ = 15 
mm, l* B = 2 mm, Db = 16 m 2 /s and T a = 160 /is. For com- 
parison, the dashed blue line shows the long-time behaviour 
predicted by diffusion theory. 

g > 1, reveals a universality of the statistics of localized 
waves. 

Our discovery of 3D transverse localization provides a 
powerful new approach for guiding future investigations 
of localization for any type of wave - not only for assessing 
whether or not the waves are localized but also for mea- 
suring the localization length. By studying three differ- 
ent fundamental aspects of Anderson localization simul- 
taneously, we have also shown the versatility of ultrasonic 
experiments, which are very well suited for undertaking 
a complete study of this phenomenon - one of the most 
fascinating in all condensed matter. Our results suggest 
that it should now be feasible to determine critical ex- 
ponents, by measuring the variation of the localization 
length with frequency near the mobility edge, as well as 
to investigate the sensitivity of the localization length to 
sample microstructure, the spatial correlations of local- 
ized waves in 3D, and the microscopic nature of localized 
states. All of these experiments, either with sound or 
light, would make unique contributions to advancing our 
understanding of Anderson localization. 
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FIG. 4: Statistical approach to localization. Comparison of the near field speckle patterns I(x,y)/ (I) for diffuse and 
localized waves observed at frequencies of 0.20 (a) and 2.4 MHz (b). In a, the speckle pattern is typical for the diffuse regime, 
while panel b reveals the narrow intense spikes that we explain in terms of Anderson localization. Note the different colour scales 
in the two figures, c, The measured probability distribution P(I) of normalized intensity I — I / (I) at 0.2 MHz (open circles) is 
close to the Rayleigh distribution (dashed blue line). The solid magenta curve is best fit to the theory of Ref.^ with g — 11.4. 
d, At 2.4 MHz, P(I) (solid symbols) shows very large departures from the Rayleigh distribution. The solid curve shows the 

theory^ with g = 0.80. At large / > 25, the data can also be described by a stretched exponential P(I) oc exp(— 1\jg~I ) with 

the same value of g (dotted curve). The large fluctuations {I 2 } = 2.74 and the large deviation from Rayleigh statistics with 
g < 1 support our conclusion that Anderson localization of sound has set in at frequencies near 2.4 MHz. 



METHODS 

The elastic network of aluminum beads was created by 
precisely controlling the flux, alloy concentration and tem- 
perature during brazing to form elastic bonds between the 
beads while preserving their spherical shape. The samples 
were made in the form of disc-shaped slabs, so as to min- 
imize edge effects. After brazing, the samples were lightly 
polished to ensure that opposite faces of the slabs were flat 
and parallel, and carefully cleaned to remove any residue of 
the brazing process. The samples were waterproofed with 
very thin plastic walls to enable pulsed ultrasonic immersion 
transducer techniques to be useo"', thereby ensuring that the 
samples remained dry when immersed in a water tank be- 
tween the generating and detecting transducers. The coher- 
ent pulse, from which v p , v g , and £ were determined, was 
then measured by scanning the position of the sample in a 
plane parallel to the sample surfa ce and averaging the trans- 
mitted field over all positional 21 . These measurements were 
made using large-diameter immersion transducers to aid in 
the spatial averaging of the transmitted field and hence in 
the extraction of the coherent component. The ratio of the 
Fourier transforms of the transmitted and input signals gave 
the amplitude transmission coefficient. 



To measure I(p,t), a quasi-point source was created by fo- 
cusing the pulsed ultrasonic beam onto a small aperture cut 
in the tip of an acoustically-opaque cone-shaped screen. The 
cone shape was chosen so that edges of the beam could be 
effectively blocked when the aperture was placed close to the 
sample being investigated, while at the same time prevent- 
ing significant stray sound being reflected back towards the 
sample from the screen. The transmitted field at a transverse 
distance p was measured using a miniature hydrophone with 
a diameter of 400 /im, which is smaller than a wavelength, al- 
lowing the transmitted field to be detected in a single speckle 
spot. The average transmitted intensity I(p,t) was deter- 
mined for selected values of p by scanning the sample posi- 
tion in a grid over a very large number of independent speckle 
spotP^. To measure the time-dependent transmission I(t) for 
an extended quasi-plane- wave source, the sample was placed 
deep in the far field of a small-diameter planar transducer, 
and the transmitted speckle pattern was scanned by moving 
the hydrophone with the sample fixed in position. For mea- 
surements of both I(p, t) and I(t), the number of speckle spots 
over which the intensity was averaged was typically 529-3025, 
and the bandwidth was set at 5% of the central frequency of 
the pulse by digitally filtering the transmitted field before de- 
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termining the dynamic intensity profiles. 

The normalized intensity I = 1/ (I) at a particular fre- 
quency was determined from the square of the magnitude 
of the Fourier transform of the transmitted field in each 
near-field speckle, normalized by the average intensity in the 
speckle pattern. The results were then binned to determine 
P (I) . To improve the statistics, the results were averaged over 
up to 100 frequencies within a 5% bandwidth where the statis- 
tics were similar — the same bandwidth as for the dynamic 
measurements of I(p,t) and I(t). In the localized regime, 
where fluctuations are largest, the statistical accuracy of the 



measured distribution was further improved by averaging over 
different samples, which were found to exhibit the same statis- 
tics. This allowed the error bars in Fig. [4ji to be determined 
from the standard error in the mean. 
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